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Abstract 

We consider a Fisher-KPP equation with density-dependent diffu- 
sion and advection, arising from a chcmotaxis-growth model. We study 
its behavior as a small parameter, related to the thickness of a diffuse 
interface, tends to zero. We analyze, for small times, the emergence of 
transition layers induced by a balance between reaction and drift ef- 
fects. Then we investigate the propagation of the layers. Convergence 
to a free-boundary limit problem is proved and a sharp estimate of the 
thickness of the layers is provided. 

Key Words: density-dependent diffusion, Fisher-KPP equation, chemo- 
taxis, drift effect, singular perturbation!^ 



1 Introduction 

In this paper we consider a Fisher-KPP equation with density-dependent 
diffusion and advection, namely 

Ut = eAiu"^) - V • (uVv^) + -u{l - u) in (0, oo) x Q 
(P^) <;^ = o ' on(0,oo)xaO 

^u{0, x) = uo{x) in , 



^ The first author is supported by the French Agence Nationale de la Recherche within 
the project IDEE (ANR-2010-0112-01). 

^AMS Subject Classifications: 35K65, 35B25, 35R35, 92D25. 
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with e > a small parameter and v^{t,x) a smooth given function. Here 
is a smooth bounded domain in {N > 2) , is the Euclidian unit normal 
vector exterior to 90 and m > 2. We are concerned with the behavior of 
the solutions u^{t,x) as e — > 0. 

Assumption 1.1 (Initial data). Throughout this paper, we make the fol- 
lowing assumptions on the initial data. 

(i) Let VIq he a nonempty open hounded set with a smooth houndary and 
such that JIq C 0. Let uq : ^ ^ he in ct^^ iiT- strictly 
positive on and such that uq{x) = for all x € d^Q. Define the 
map uo : — > K by 



Uo{x) :-- 



uo{x) if X G flo 

if X e Q\Ti^ . 



(a) Qq is convex. 

(Hi) there exists 5 > such that, if n denotes the Euclidian unit normal 
vector exterior to the "initial interface" Tq := dil.(), then 



duo, . 



> 5 for all y G To 



(1.1) 



Remark 1.2. Note that the comparison principle allows to relax the regular- 
ity assumption on uq- See [2] for details. 

Assumption 1.3 (Structure of v^). We assume that 

v''{t,x)=v{t,x)+£vl{t,x), (1.2) 

with V and vf smooth functions on [0, oo) x 0. We assume that, for all 
T > 0, there exists C > such that, for all e > small enough, it holds that 
Ikf llci'2([o T]xn) — ^- Finally we assume 

— = on {0,oo)x on. (1.3) 

Remark 1.4. In the sequel we smoothly extend v{t,x) in time and space on 
the whole of R X R^, as well as vl{t,x) in space on [0, cxd) x R^. Moreover 
since we are investigating local in time phenomena we will assume in the 
sequel, without loss of generality, that the extensions v{t,x) and uf(t,x) 
vanish outside of a large time-space ball. 

Problem (-P^) is a simpler version of a chemotaxis-growth system with 
logistic nonlinearity, where v^{t,x) is not a given function but is coupled 
to u either through the parabolic equation evt = Av + u — ^yv ot through 
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the elliptic equation = Av -\- u — ^v, supplemented with the Neumann 
boundary condition (jl.3p (see e.g. [21j). Note that, in the case of linear dif- 
fusion (corresponding to m = 1) and a bistable nonlinearity, the asymptotic 
behavior of the corresponding system as e ^ has been studied using the 
Green's function associated to the homogeneous Neumann boundary value 
problem on for the operator — A + 7 (see [8] and [1]) . 

Motivation and biological background. Before describing our results, 
let us briefly comment about the relevance of {P'^) in population dynamics 
models. The evolution equation in Problem (-P^) combines logistic growth, 
chemotaxis and degenerate diffusion. We recall below how these terms ap- 
pear in mathematical models that attempt to capture remarkable biological 
features. 

Reaction diffusion equations with a logistic nonlinearity were introduced 
in the pioneering works [12], [18]. The simplest equation reads 

ut = Au + n(l — u) , 

and has been widely used to model phenomena arising in population genetics 
[12j or in biological invasions [22]. Its main mathematical property is to 
sustain travelling wave solutions with a semi-infinite interval of admissible 
wave speeds, with the minimal one having a crucial biological interpretation. 

Chemotaxis, i.e. the tendency of biological individuals to direct their 
movements according to certain chemicals in their environment, is induced 
in (P^) by the advection term —V • (uVf*^): the population, whose density 
is u{t,x), has an oriented motion in the direction of a positive gradient 
of the chemotactic substance, whose concentration is v'^{t,x). The first 
PDE model to describe such movements was proposed in [T7] and involves 
linear diffusion for u and a parabolic equation coupling v to u. The Keller- 
Segel model has received considerable attention in mathematical literature, 
particularly focusing on the finite-time blow-up of solutions (see [16] for a 
recent review). This provides a mathematical tool to analyze aggregation 
phenomena as observed in bacteria colonies. Chemotaxis systems involving 
linear diffusion and a growth term, either logistic or bistable, have later been 
considered in, e.g., [21], [8], [1] and [21] . 

Variants of the Fisher-KPP equation involving a degenerate diffusion 
have been proposed in order to take into account population density pres- 
sure. Actually one can introduce density-dependent birth or death rates as 
an attempt to control the size of a population. Nevertheless as shown in |13j . 
the introduction of a nonlinearity into the dispersal behavior of a species, 
which behaves in an otherwise linear way, may lead, in an inhomogeneous 
environment, to a similar regulatory effect. Moreover this assumption is 
consistent with ecological observations as reported for instance in [S], where 
it is shown that arctic ground squirrels migrate from densely populated areas 
into sparsely populated areas, even when the latter is less favorable (due to 
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reduced availability of burrow sites or exposure to intensive predation). For 
such species, migration to avoid crowding, rather than random motion, is 
the primary cause of dispersal. To describe such movements, the authors in 
|22j and [13] use the directed motion model where individuals can only stay 
put or move down the population gradient; this model yields the degenerate 
equation 

ut = A{u^) + G{x)u , (1.4) 

in which the population regulates its size below the carrying capacity set 
by the supply of nutrients. Later in |14j a larger class of equations with 
degenerate diffusion and nonlinear reaction was considered, namely 

ut = A{u'^) + f{u), m>2. (1.5) 

Note that in the absence of f{u), equation p.Sp reduces to the so-called 
porous medium equation 

ut = Aiu"'), (1.6) 

which has been extensively investigated in the literature. We refer to the 
book [23] and the references therein. The main feature of this equation is 
that it is degenerate at the points where u = 0. As a consequence, a loss of 
regularity of solutions occurs and disturbances propagate with finite speed, 
a property which has a relevant interpretation in a biological context (see 
for instance [6]). 

Formal asymptotic analysis. Problem (-P^) possesses a unique solution 
u^{t,x) in a sense that is explained in Section [3l As e — )• 0, the qualita- 
tive behavior of this solution is the following. In the very early stage, the 
nonlinear diffusion term eA(n'") is negligible compared with the drift term 
— Vu • Vu^ and the reaction term e~^u{l — u). Hence, in some sense, the 
equation is well approximated by a coupling between the transport equation 
ut + Vu ■ S/v^ = and the ordinary differential equation ut = e"^u{l — u). 
Therefore, as suggested by the analysis in [2], u*^ quickly approaches the 
values or 1, and an interface is formed between the regions {u^ ~ 0} and 
{n^ ~ 1} [emergence of the layers). Note that, in this very early stage, 
the balance of the transport equation and the ordinary differential equation 
will generate an interface not exactly around Fq but in a slightly drifted 
place. Once such an interface is developed, the diffusion term becomes large 
near the interface, and comes to balance with the drift and the reaction 
terms so that the interface starts to propagate, on a much slower time scale 
[propagation of the front) . 

Our goal in this paper is to provide a rigorous analysis that supports 
this formal approach and makes it more precise. To study the interfacial 
behavior, we consider the asymptotic limit of (-P^) as e 0. Then the limit 
solution will be a step function u[t, x) taking the value 1 on one side of a 
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moving interface, and on the other side. We show that this sharp interface, 
which we will denote by Tt, obeys the law of motion 



where Vn is the normal velocity of Tt in the exterior direction, c* the minimal 
speed of travelling waves solutions of a related degenerate one-dimensional 
problem (see Section [5] for details) and n the outward normal vector on Tt- 

Plan. The organization of this paper is as follows. We present our results in 
Section [5J In Section [21 we briefly recall known results concerning the well- 
posedness of Problem (P"^); in particular, it admits a comparison principle 
so that the sub- and super-solutions method can be used to investigate 
the behavior of the solutions u^. In Section U we prove a generation of 
interface property for Problem (P^). In Section[5]we investigate the motion 
of interface. Finally, we prove our main result in Section [6l 

2 Results and comments 

The question of the convergence of Problem (P^) to (P^) has been addressed 
in [TT]. However, the author considers only a very restricted class of initial 
data, namely those having a specific profile with well-developed transition 
layers. In other words the generation of interface from arbitrary initial data 
is not studied. In the present paper we study both the emergence and the 
propagation of interface. Moreover we prove a sharp 0{e) estimate of the 
thickness of the transition layers of the solutions . 

The authors in [15] prove the convergence of the solutions of (P^) with 
arbitrary initial data with convex compact support to solutions of (P*^), 
when there is no advection (i.e. = 0). They provide an 0{e\ In e\) estimate 
of the thickness of the transition layers. Therefore, even in the particular 
case = 0, our 0(e) estimate was not known. 

As mentioned in the introduction, the drift term and the reaction term 
in (P^) are of the same magnitude for small times. Therefore the emergence 
of the layers, initiated by the ODE ut = s~^u{l — u), will occur in the 
neighborhood of a slightly drifted initial interface Fq To analyze such a 
phenomenon we shall use the Lagrangian coordinates. Recall that we have 
smoothly extended v{t, x) in time-space on the whole of M x M^, with v = 
outside of a large time-space bah. Then, for {to,xo) e M x M^, we denote 
by (p[to,xo) the solution, defined on M, of the Cauchy problem 





(2.1) 
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We denote by $ the associated flow defined on M x M x M'^, that is 

«>(ti,t2,X3) :=^(t,,.3)(ti). (2.2) 
Recall that Tq = dilo = 9(Supp uq) is the initial interface. From i = to 

:= e\ lne| (generation time) , (2.3) 



we let each point on Tq evolve with the law ()2.ip and then define a drifted 
initial interface Tq " by 

rf"f*:={$(t%0,x): xGTo}. (2.4) 
Next we consider the free boundary problem 



(-f'e'drift) ' 



dv 



Vn = c* + — on 
on 

pE, drift I pE, drift 

\t=0~^0 



£, drift 



Well-posedness of (P^) and of (^'g drift)* Using the level set formulation 
(see, e.g., [5]), the motion law in Problem {P^) can be rewritten as a first 
order Hamilton-Jacobi equation with a convex Hamiltonian. This approach, 
combined with the results in [19] , has been used in [11] in order to prove the 
following. 

Theorem 2.1 ([11], Well-posedness of (-P°)). Let CC Q be a smooth 
subdomain of and let Tq = d^o be the given smooth initial interface. 
Then there exists T"^'^^{Tq) > such that Problem (P^) has a unique smooth 
solution on [0,T] for any < T < T'""^'(ro). More precisely, there exists a 
family of smooth subdomains {^t)tG{o,T] ''JJ'i-th CC O such that, denoting 
Tt = dQt, r := Uo<i<T({*} ^ ^t) ^-^ unique solution to Problem {P^) on 
[0,T]. 

Moreover, r"^"^(ro) depends smoothly on Tq. Therefore we can choose 
eo > small enough and T > such that 

0<r< inf T"*'^^(r^''^"^*) , (2.5) 

0<£<£o 

which guarantees the existence of a unique smooth solution on [0, T] to both 
(PO) and (if drift) for any < e < eo- We denote by T^'^^^f* = Uo<t<r({*} ^ 
-p£,drift-j ^j^g smooth solution to (if drift) and by Qf'^'^^^^ the region enclosed by 
-pe,drift^ In the sequel we work on [0, T], with T satisfying (j2.5|) . and define 

Qt := (o,r) X n. 

Our main result. Theorem 1 2. 2 1 contains generation, motion and thickness 
of the transition layers properties. It asserts that: given an initial data uq, 
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the solution quickly (at time = e| lne|) becomes close to 1 or 0, except 
in a small neighborhood of the drifted interface Fji , creating a steep 
transition layer around p^j'^^''^* [generation of interface). The theorem then 
states that the solution remains close to the step function associated with 
(-^edrift) time interval [motion of interface); in other words, 

the motion of the transition layer is well approximated by the limit interface 
equation (-Pg drift)- Moreover, the estimate (|2.6p in Theorem 12.21 implies that, 
once a transition layer is formed, its thickness remains within order 0[e) for 
the rest of time. 

Theorem 2.2 (Generation, motion and thickness of the layers). Let r] G 
(0, 1/2) be arbitrary. Then, there exists C > such that, for all e > small 
enough and all 

f = ellnel <t<T, 



we have 



f [0,1 + 77] if xeMceiT] 



u'^[t, x) G < 



£,drift^ 



[1-7?, 1 + 7?] Z/ XGnf'"^'\Mce{rf''"^) (2.6) 



Uo} z/ xG[n\nt'''"^')\Mce[rf''"^) 

with Mr [Tf'^"-^'') := {x : dist(3;, Fj''^"-^*) < r} the tubular r -neighborhood of 



Y^, drift 

Note that ([221) shows that, for any < a < 1, for aU t" < t <T, the 
a-level set 

Lf(a) := {x : u^[t,x) = a} 

lives in a tubular 0[e) neighborhood of the interface y^''^^^^^. In other words, 
we provide a new 0[£) estimate of the thickness of the transition layers of 
the solutions u'^ . Concerning the localization of the level sets L^[a), it is 
made with respect to a slightly drifted free boundary problem (-P^jj-ift)- 
Nevertheless, since the solution of (-P(? drift) [0)^] is continuous w.r.t. the 
initial hypersurface Fq'^"^*, we recover, as e — > 0, the original free boundary 
problem [P^) and obtain the expected result. More precisely, let us define 
the step function u[t, x) by 

u[t,x):=\ . \ fortG(0,r]. (2.7) 

As a consequence of Theorem 1 2. 2 1 we obtain the following convergence result 
which shows that u is the sharp interface limit of as e — > 0. 

Corollary 2.3 (Convergence). As e ^ ^, u'^ converges to u, defined in 
(|2.7p . everywhere in Uo<t<T({*} ^ ^t) Uo<£<r ({*} ^ {^\ ^t)) ■ 
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3 Comparison principle, well-posedness for (P^) 



Since the diffusion term degenerates when n = a loss of regularity of 
solutions occurs. We define below a notion of weak solution for Problem 
(P^), which is very similar to the one proposed in [3] for the one dimensional 
problem with homogeneous Dirichlet boundary conditions. Concerning the 
initial data, we suppose here that uq € L°°{Q) and uq > a.e. Note that in 
this section, and only in this section, we assume, for ease of notation, that 
e = 1 and that = v; we then denote the associated Problem (P^) by (P). 
In the sequel f{u) = u{l — u). 

Definition 3.1. A function u : [0, oo) — )• L^{Q) is a solution of Problem 
(P) if for all T > 0, 

(i) u £C {[0,oo);L^n)) nL^^iQr) 

dip 

(a) for all ip € C'^{Qt) such that (p > and = on dQ, it holds that 

/ u{T)ip{T) -11 {uipt + u'^Aif + uVv ■ Vif) 
Jn J Jqt 

= I uom+ I [ fiu)^- (3.1) 
Jn J Jqt 

A sub-solution (a super-solution) of Problem (P) is a function satisfying (i) 
and {ii) with equality replaced by < (respectively >). 

Theorem 3.2 (Existence and comparison principle). Let T > be arbitrary. 
The following properties hold. 

(i) Let u~ (w^ ) be a sub-solution (respectively a super- solution) with ini- 
tial data Uq (respectively Uq). 

If Uq < Uq a.e. then u~ < in Qt; 

(ii) Problem (P) has a unique solution u on [0, oo) and 

< ti < max(l, ||tio||ioo(Q)) inQx; (3-2) 

(Hi) u G C(Q^). 

Since (jl.Sp holds, the proof of Theorem 13.21 is standard and follows the 
same steps of that of [3l Theorem 5] . The continuity of u follows from [10] . 

The following lemma proved in [15], will be very useful when constructing 
smooth sub- and super-solutions in later sections. 
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Lemma 3.3 (Being sub- and super-solutions). Let u be a continuous non- 
negative function in Qt- Define 

Qf'VP ,= {xen: u{t, x) > 0} rfPP := 9^!^^^ , 

for all t G [0,r]. Suppose the family T'^^pp := Uo<t<T ({*} x ^T^^) suffi- 
ciently smooth and let j/^"^^ be the outward normal vector on Y^^^^. Suppose 
moreover that 

(i) V(u'") is continuous in Qt 

(ii) C''[u]:=ut-A{u'^) + V-{uVv)-f{u) = in {{t,x) eQ^ : u{t,x) > 
0} 

(in) ■?S = on 517^^^ for all te[0,T]. 

t 

Then u is a solution of Problem (P) . Similarly a function satisfying (i) and 
(a) — (Hi) with equality replaced by < (>) is a sub-solution (respectively a 
super-solution) of Problem (P). 



4 Emergence of the transition layers 

In this section, we investigate the generation of interface which occurs very 
quickly around F^e'^'^'^*. We prove that, given a virtually arbitrary initial 
datum uq, the solution of (P^) quickly becomes close to 1 or in most 
part of $7. More precisely — recalling that ^{ti,t2, x^), defined in ()2.2p . 
denotes the flow associated with the Cauchy problem (j2.ip — the following 
holds. 

Theorem 4.1 (Emergence of the layers) . Letrj G (0,1/2) be arbitrary. Then 
there exists Mq > such that, for all e > small enough, the following holds 
with = e| lne|. 

(i) for all X £ ^l, we have that 

< M^(^^x) < 1 + (4.1) 

(a) for all X G Q, we have that 

if uo{^{0,t^,x)) > Moe then ^^^(^^ x) > 1 - ry ; (4.2) 

(Hi) for all X G Q, we have 

if dist($(0,^^x),^^o) > Moe then n^(^^x) = 0, (4.3) 

where we recall that Qq = {x : uq{x) > 0} (see Assumption \l.l\} . 

In order to prove the above theorem, we shall construct sub- and super- 
solutions. As mentioned before, in this very early stage, we have to take 
into account both the reaction and the drift terms. We start with some 
preparations. 
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4.1 A related ODE and the flow $ 

An ODE. The solution of the problem without diffusion nor advection, 
namely ut = f{u) supplemented with the condition u{0,x) = uo{x), is 
written in the form u{t,x) = Y uq{x)^, where y(r, ^) denotes the solu- 
tion of the ordinary differential equation Yt-{t, ^) = /(y(r, ^)) supplemented 
with the initial condition Y{0, ^) = ^. Nevertheless, in order to take care of 
the term —uAv^, we need a slight modification of /. 

Let / be the smooth odd function that coincides with f{u) = n(l — u) 
on [0, oo): / has exactly three zeros — 1 < < 1 and 

/'(-1) = -1<0, /'(0) = 1>0, /'(1) = -1<0, (4.4) 

i.e. / is of the bistable type. Next, we define 

Mu) ■.= f{u) + 5. 

For \5\ small enough, this function is still of the bistable type: if 5o is small 
enough, then for any 6 G {—Sq, 6o), fs has exactly three zeros a-{6) < a{6) < 
a+{6) and there exists a positive constant C such that 

\a^{6) + 1| + \a{6)\ + |q+(5) - 1| < C\6\ , (4.5) 

H6)-1\<C\5\, (4.6) 
where //(J) is the slope of fs at the unstable zero, namely 

H{6) := /Ka(5)) = f'{a{6)) . (4.7) 

Now for each 6 € {—6o,6o), we define Y{t,S, ; (^) as the solution of the 
ordinary differential equation 

Yr{T,(;5) =MY{t,C,S)) forr>0 

where ^ varies in (—Co, Co), with 

Co := ||^^o||l-'{C) + 1 • (4.9) 
We claim that Y(t,S^; 5) has the following properties. 

Lemma 4.2 (Behavior of Y). There exist positive constants 6q and C such 
that the following holds for all (r, ^ ; 5) € (0, oo) x (—Co, Co) x {—5q, 5q). 

(i) If i>a{5) then Y{T,i-5) > a{5) 
If i<a{5) then Y{T,i]5) < a{5) 

(n) \Y{T,i-5)\<Co 
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(Hi) <Y^{t,^;6) < Ce^(^)" 



(iv) 



< C{e^'^^> - 1). 



Properties (i) and (ii) are direct consequences of the bistable profile of 
fs- For proofs of {in) and (iv) we refer to [1]. 

The flow <1>. Let us briefly recall well known facts concerning the flow 
^(ti,t2, x^). By definition we have 



5$ 



it,to,xo) = Wv{t,<^{t,to,xo)) . 



(4.10) 



Next, note that, by uniqueness, 

<^{t,to,'^{to,t,xo)) = xo, 

for all {t,tQ,XQ) € M X M X M.^ . Differentiating this identity with respect to 
to, we get 

— {t,to,x) + D3<^{t,to,x) — {to,t,Xo) = OjjiV , 

where x := ^{tQ,t,xo) and where D3^{ti,t2, x^) denotes the Jacobian ma- 
trix of $ w.r.t. the third variable. Hence, using ()4.10p we infer that 



— (t, to,x) + Ds^it, to,x)Vv{to,x) = OjRiv , 
which is of crucial importance for our analysis. 



(4.11) 



4.2 Proof of (KB and (K2 



We use the notation z'^ = max(2;, 0). Our sub- and super-solutions are given 
by 



{t, x) : 
or equivalent ly by 

wf{t,<^{t,0,x)) : 



Y ( ^, no($(0, t, x)) ± e2^*(eA'(±^^)*/^ - 1) ; ±eM 



(4.12) 



Y (^^, uo{x) ± e2^*(eM(±^^^)*A _ i) ; ±eM^ 



(4.13) 

Here Y{t,S,;6) is the solution of (j4.8p . fi{6) the slope defined in (j4.7p . 
^{ti,t2,X3) the flow defined in ()2.2p and M is chosen such that, for all 
e > small enough, M > Co||^'y^||L°o(Qy)i with Cq defined by ()4.9p . 
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Lemma 4.3 (Sub- and super-solutions for small times). There exists C* > 
such that, for all e > small enough, {w~,wf) is a pair of sub- and super- 
solutions for Problem {P'^), in the domain [0, i^] x Q,. 

Before proving the lemma, we remark that w~{Q, x) = w^(0, x) = uq{x). 
Consequently, by the comparison principle, we have 

w-{t,x) <u''{t,x) <wt{t,x) for all G [0,t^] X 0. (4.14) 

Proof. In order to prove that (w~, wf) is a pair of sub- and super-solutions 
for Problem {P'^) — if C* is appropriately chosen — we check the sufficient 
conditions stated in Lemma 13.31 

On the one hand, concerning the sub-solution to", for (t,x) such that 
X G VlI^^^[wJ] := {x : w~{t,x) > 0} we have, at point {t,x), 

V{w~r = (my^-iyg)^^, no($(0,t,x)) -e2c7*(e'^(-^^)*/" - 1); 

-eM)v^(«o($(0,t,2;)). 

If {t,x) — > {to,XQ) such that xq G r*^^^[^y~] := d^ll^^[w~] then the equality 
above implies 

lim V{w-r{t,x) = O^N , 

and conditions (i) and (iii) of Lemma 13.31 are checked for the sub-solution. 

On the other hand, concerning the super-solution , note that ^ := 
uo{^{0, t, x)) -\- e2(^*(^gM(£A^)t/e _ -^-^ jg positive. Therefore the cubic profile of 
fs shows that, for t > 0, 

r'^''pp[w+] := = dn . 

Recall that uq = in a neighborhood V of dO,; if x is sufficiently close to 
dQ, ^{0,t,x) lives in V for all t G [0,t^] (with e > sufficiently small). 
Therefore ()4.12p shows that is independent on x near and condition 
(Hi) of Lemma 13.31 for the super-solution is checked (and condition (i) is 
obviously checked). 

Then it remains to prove that 

£^[w~] := iw-)t - sAiw^r + V • {w^Vvn - -Jiw^) < , 

in {{t,x) G [0,^] X n : wj{t,x) > 0} and that C^[ivf] > in {{t,x) G 
[0,t^] X 0}. We will only prove the latter inequality since the proof of the 
former is similar. 
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We compute 

dtwt = U'r + ^^ [^xo(a>(0,t,x))] + e/i(eM)C*e^(^*W^y5 

Vu;+ = V,[no(«5(0,t,x))]yj 
V[{wtr] = V,. [uo(a>(0,t,x))] (y™)g 

A [{wtT] = I V,. [no(^>(0, t, x))] P(y™)55 + A,. [no(<^(0, x))] (r-)^ , 
where the function Y and its derivatives are taken at the point 

(r, i-6):= {t/e, uo{m t, x)) + e'C^(e'^(-^We _ . . 
Note that 

d 

— [uomO,t,x))] = Vuoi^{Q,t,x)) ■ Q^iO,t,x) , 

and that 

[no($(0,t,x))] = p3$(0,t,x))^Vno($(0,t,x)) , 

with (Z)3<I>(ti, t2> a^s))^ the transpose of the Jacobian matrix of <I> w.r.t. the 
third variable. 

Therefore, using f{wf) = f{wf) = fehiiY) — eM and the equation 
Yt = feM{y), we infer that 



C'[wt]=Ei+E2 + eY^Es, 



where 
El 



M + YAv'' 



E2 = VuoiHO, t, x)) ■ (||(0, t, x) + D3$(0, t, x)Vvit, x)) Y^ 

E3 = C*fi{eM) e/^(^^We + {Ds^iO, t, x))^Vuo($(0, t, x)) ■ Vvl{t, : 



-A^ [uoi^iO,t,x))] 



|V, [no($(0,t,x))]|2iI^ 



We note that, for e > sufficiently small, 6 = eM € {—So, 60) and that, 
in the range < t < = e| lne|, 

e = no(<I>(0,t,x)) +e2c7^(eMMW. _ 1) ^ (-Co,Co) , 

so that estimates of Lemma 14.21 on Y{t,S,;6) will apply. 

Since we have chosen M > Co||At''^||^oo(Q^), Ei > holds. Moreover, 
(j4.1ip implies E2 = 0. In the sequel we denote by C various positive con- 



stants which may change from place to p 



Lemma IT2] (u) — (iv) we see that 



ace but do not depend on e. From 
= Imy^-^l < C and that 



Y, 



< m(m - i)|y™-2y.| + my™-^ 



Y 



55 



Y, 



<C + C{e 
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since m > 2. Hence 

E3 > (CV(eM) - C7)e'^(^A./)Ve _ ^ _ 

Since iJ,{eM) 1 as e ^> 0, by choosing C* » C we see that E-^ > for all 
e > small enough. 

Recalling that > 0, we get > and the lemma is proved. □ 

We are now in the position to prove ()4.ip and ()4.2p . 

Proof. Let r] € (0,1/2) be arbitrary. Then [H Lemma 3.11] provides a 
constant Cy > such that, for all e > small enough, for all ^ G (— Cq, Cq), 

y(|lne|,^;±eM) < 1 + 7?; (4.15) 

if C>Cye then y(| lne|, ^ ; ±eM) > 1 - r/ . (4.16) 
By setting t = = e| In e| in (|4.14p . we obtain 

Y [\ lne|, no($(0, f,x)) - CVe^^"^^)! '"^^I + ; -em) ^ < u^f, x) 
<y(|lne|,no($(0,^^x)) + CVe'^(^*^)|l°^l-CV;eM)^ . (4.17) 

Therefore, the assertion ()4.ip of Theorem 14.11 is a direct consequence of 
KT7\i and (015]). Next we prove (jM]). Note that in view of (gS]), we 
have se^^~'^^^^^^^^ — >• 1 as e ^ 0. Therefore, for e > small enough (since 

n>o), 

«^(^^2;) >y^|lne|,no($(0,^^x))-^C*e + CV;-eM^ . (4.18) 

Choose Mo > so that MqE - ^C* e + C* e"^ > max(Cye, a(-eM)), with Cy 
as in (|4.16p . Then, for any x G i7 such that tio($(0, t^, x)) > MqE, we have 

«o(^(0, x)) - ^C*e + CV > Cye . 

Combining this, ()4.18p and ()4.16p . we see that 

u%f,x) > 1 - ??. 

This completes the proof of ()4.2p . □ 
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4.3 Proof of 

Let us recall that a finite speed of propagation property, as is ()4.3p . is proved 
in [15]: the authors construct a super-solution using a related travelling wave 
U of minimal speed, and they obtain an ©(el lne|) estimate of the thickness 
of the transition layers. We borrow some ideas from this paper but, in order 
to obtain the improved 0{e) estimate, we again use the solution Y of the 
ordinary differential equation ()4.8p . 

Let be the solution of the Cauchy problem (recall that v'^{t,x) has 
been extended on [0, oo) x in Remark II. 4|) 

f = eA(z"^) - V • (zVv^) + -f{z) in (0, oo) x 
[z{0,x) = uo{x) inR^. 

Lemma 4.4 (Super-solutions for (Q^) for small times). Choose K > 1 and 
> appropriately. For all xq € OQq = 9Supp uq, denote by uq the unit 
outward normal vector to dilo at xq. For t > 0, x € R", define the function 

z^{t,x) := 



K 



Y (^-, -($(0,t,x) - xo) • no + e'a(e^(^*-^)*/^ - l);eM^ 



Here Y{t,^]5) is the solution of ()4.8p . /u(J) the slope defined in ()4.7p . 
$(ti,t2ia;3) the flow defined in (j2.2p and M is chosen such that, for e > 
small enough, M > Co|| Af^||/^oo(Qy) . Then, for all e > small enough, 

Uo{x) < 4(0^ ^) for all x G R^ , (4.19) 

and 

C^[zt] := izt)t - sAiztr + V • {ztVv^ - ^/(^+) > , (4.20) 



in the domain [0, t^] x 



Proof. Recall that f^o is convex. Therefore, in view of p.ip . we can choose 
K > 1 sufficiently large so that, for all xq G OQq and all x £ ^Iq, 

uoix) < -K{x - Xq) ■ no . (4.21) 

We prove (j4.19p . If <&(0,0, x) = x ^ ^Iq this is obvious since uo{x) = 0. 
Let us now assume that <1>(0, 0,x) = x £ Qq. Since f^o is convex, it lies on 
one side of the tangent hyperplane at xq so that (x — xq) • tiq < 0. Recall 
that Y{0,^;6) = ^ so that z+(0,x) = -K{x - xq) ■ hq and follows 
from dlT^ . 
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We now prove ()4.20p . As in the proof of Lemma 14.31 straightforward 
computations combined with (|4.8p and (j4.1ip yield 

eC'[zt] = Kf{Y) - f{KY) + £K{M + Y /\v^) 
+e^KY^\c^lx{eM) e'^("^'^)*/^ - D3$(0, t, x)no ■ Vv{{t, x) 

+A,[«J>(0,t,x) . no]K™-i^^ - |V,[<^.(0,t,x) . n^]\^K^-^^^] . 

Note that Kf{Y) - f{KY) = K{K - 1)Y'^ > 0. Then, by using similar 
arguments to those in the proof of Lemma 14.31 see that > 0, if 

Cj, > is sufficiently large. □ 

We now prove (|4.3p . 

Proof. We shall first prove that property (|4.3p holds for the solution of the 
Cauchy Problem (Q^)- Recall that a{S) is the unstable zero of fs = f + 6 
so that a{eM) < 0. Moreover, in view of (j4.5p and (j4.6p . we can choose 
Mq > large enough so that, for e > small enough, 

-Moe + aee(^("*-^)-i)l '"^"1 - C^e^ < a{eM) . 

For X € such that dist(<I>(0, t^, x), f^o) ^ Afoe, we choose xq G d^Q such 
that dist($(0, t^, x), Oq) = |1^(0, ,x) — xq\\ and define zf as in Lemma 
It follows from Lemma 14.41 and the comparison principle that, for all e > 
small enough, ah {t,x) e [0, i^] x M^, 

Q < z^{t,x) < zf{t,x) . (4.22) 

Since, for t = = e\\u. e\, 

- (a>(0, ^^ x) - xo) ■ no + e2c,(eMeM)tV^ _ i) 
= -\\<^iO,f,x) - xoW + c^ee(^("*^)-i)l''^^l - 

< a{£M) , 
it follows from Lemma 14.21 (i) that 



Y ( |, -($(0, t", rr) - xo) • no + e2c^(eA'M^)*V^ - 1) ; eM ) < a(eM) < , 



and therefore zj{t^ x) = 0, which in turn implies x) = 0. Hence ()4.3p 
holds for the solution of (Q^)- 

Now, a straightforward modification of |15l Corollary 4.1] shows that 
there exists T > such that, for all e > small enough, 

x) = x) , 

for all (t,x) e (0,r) X VL. This proves ([O]) for u" the solution of (P^). □ 
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5 The propagating front 



The goal of this section is to construct efficient sub- and super-solutions that 
control n"^ during the latter time range, when the motion of interface occurs. 
We begin with some preparations. 

5.1 Materials 

In the linear diffusion case (m = 1), it is well-known that the equation 
ut = An-|-it(l — u) admits travelling wave solutions with some semi-infinite 
interval of admissible wave speed. The same property holds for the nonlinear 
diffusion case, namely equation ut = A(ti™)-|-ti(l — n), m > 1. Nevertheless, 
it turns out that the travelling wave with minimal speed c* > is both 
compactly supported from one side and sharp. In the following, U denotes 
the unique solution of 

' + c*U'{z) + U{z){l - Uiz)) = for all z G M 

= ^ (5,1) 

U{z) > for all z < 

U{z) =0 for all z > . 

Lemma 5.1 (Behavior of U). For all z € (— oo,0) we have U'{z) < 0. The 
travelling wave U is smooth outside and 



U'{0) 



= if 1 <m <2 

€ (-00,0) if m = 2 
= —oo if m > 2 . 



Moreover, there exist C > and /3 > such that the following properties 
hold. 

|(C/"^)'(z)| < CU{z) for all z£R, (5.2) 

< 1 - Uiz) < Ce^^l^l for all z < , (5.3) 

\zU'{z)\ < CU{z) for all z < -I . (5.4) 

For more details and proofs we refer the reader to [1], [7|, [15], as well 
as to Ho] for related results. 



Another ingredient is a "cut-off signed distance function" d^{t,x) which 
is defined as follows. Let d"^ = d'^''^^^^^ be the signed distance function to 
-p£,drift^ the smooth solution of the free boundary problem (-P(?drift)' namely 



~, , I -dist(x,i/ ) tor X G SL' 
^ef^ £^ ._ ) . — (5.5) 

dist(x,Fj' " ) for X € \ Oj' " 
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where dist(x,r^' is the distance from x to the hypersurface T^' We 
remark that = on r^'^^i'ift and that | Vd^ | = 1 in a neighborhood of r^''^"^*: 
there exists (Iq > such that, for all e > small enough, \S/d^{t,x)\ = 1 
if \d^{t,x)\ < 2dQ. By reducing^do if necessary we can assume that is 
smooth in {{t,x) G [0,T] x H : \d%t,x)\ < 3do} . 

Next, let be a smooth increasing function on M such that 



as) 



s if |s| < do 

-2do if s < -2do 
2do if s > 2do . 



We then define the cut-off signed distance function = d'^''^'^^^^ by 

d'{t,x) := C {d'{t,x)^ . (5.6) 

Note that 

if \d''{t,x)\ < do then \Vd%t,x)\ = l, (5.7) 
and that the equation of motion (-P^P^ji-ift) recast as 

(4 + c* + Vd' ■ Vv) {t, x) = on F^''^"^* = {x e n : d'{t, x) = 0} . 

(5.8) 

Moreover, there exists a constant D > such that, for all e > small 
enough, 

\Vd^t,x)\ + \Ad''{t,x)\ < D for all (t,x) G Q^. (5.9) 

Finally, in view of (|5.7|) and (|5.8|) . the mean value theorem provides a con- 
stant > such that, for all e > small enough, 

\dt + c*\Vd^\'^ + Vd^ ■Vv\{t,x) <N\d%t,x)\ for ah {t,x) gQ^. 

(5.10) 

5.2 Sub- and super-solutions 

We define 

ut(t..):^(l±m)u('^i^^^^^^^]. (5^11) 



where 

p{t) ■= -e-*/= + e^^ + K 

q{t) ■= (j(e"*/^ + eLe^*) , 

and where C/ and were defined in subsection l5.1l In the following lemma, 
nf'P^lufl Tf'^^lufl V^PPluf] and i/f^^in^] are defined as in Lemma [331 
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Lemma 5.2 (Sub- and super-solutions for the propagating front). Choose 
o" > small enough so that 

c*{m-l)D\l + 2ar-^a (5.12) 

where D is the constant that appears in (|5.9p . Choose K > 1. Then, if 
L > is large enough, we have, for e > small enough, 

Ciu"] <0 < C'[uf] in{0,T)xn, (5.13) 

. U -1 = Q dnr'K] for all te[0,T], (5.14) 

•^^"^^ =0 on for all te[0,T]. (5.15) 



Proof Properties (|5Ji]l and (l5J5]l follow from ([/™)'(0) = (see (fOjl ). 
We prove below that >C^[u^] > 0, the proof of >C^[n~] < following the 
same lines. Note that we only need to consider the case d^{t,x) < £p{t) 
since, if (P{t,x) > ep{t) then M+(f,x) = 0. Straightforward computations 
and equation (|5.ip yield 

e{u+)t = eq'U + (1 + q)idf - ep')U' 
eVut = {l + q)U'Vd^ 

e2A(n+)"* = (1 + g)™|Vd"|2(-c*C/' - U{1 - U)) + e(l + g)'"Ad=([/'")' , 
where arguments are omitted. Thus we get 

eC'[ut]=Ei+E2 + E^, 

where 

El = U'{1 + q) [dl - ep' + c*(l + qT-^\Vd^^ + Vd^ • Vu=] =: U\\ + q)E\ 
£^2 = f/ {-(1 + g) + (1 + g)'"|Vd^|2 + [/ [(1 + qf - (1 + g)™|V(i^|2] + eg'} 
^3 = -e(l ^ q)"^Ad%U"'y + e(l q)Av'U . 

In the sequel we define a{t) := 1 -|- q{t) and denote by Q various positive 
constants which do not depend on e. 

Since ||At;'^||j;^oo(Qj,) is uniformly bounded w.r.t. e > (see Assumption 
fOD . we deduce from and that jE^sj < eCsia"' + a)U so that 

^2 + ^3 > - a + a"" + [/ (a2 - a"") - eCsa"" - eCsa 

+ (|Vd"|2-l)a™(l-[/) + eg'}. (5.16) 
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We claim that, for e > small enough, 



Kivrf 



l){l-U)\<eC2. 



(5.17) 



Indeed, if —do < d^{t,x) < £p{t), it follows from ()5.7p that, for e > small 
enough, |Vd^(t, x)| = 1. Next, if d^{t,x) < —do, ()5.3p implies that 



and (|5.17p holds for e > small enough. Therefore ()5.16p and ()5.17p imply 



E2 + E3>U{-a + a"' + U{a^- a"") - e{C2 + C3)a'" - eC^a + eq'] . 



Next, since 

E{ = dl + c*|Vd^|2 + V(f ■ Vv - ep' + c*(a™~^ - l)\Vd^\^ + eVd^ • Vv{ , 
using (|5.1Up , (|5.9p and Assumption 11.31 we see that 

E{ < N\d'{t,x)\-ep'{t) + c*{a'''-'^ -l)\Vd'f + £Vd' -Vvf 

< N\d^t,x) - £p{t)\ + £{Np{t) - p'{t)) + c*(a™~^ - 1)1)2 ^ ^^j^ 

< N\d'it,x)-£pit)\+£iNpit)-p'it)) 
+c*{m - 1)1)2(1 + 2a)™-2g(t) + eCD . 

The last inequality above comes from the fact that, for e > small enough, 
we have < q(t) < a{l + eLe^'^) < 2a, which in turn implies that 



In the following, we distinguish two cases. 

First, assume that < d^{t,x) < £p{t) so that, for e > small enough. 



El < e{2Np{t) - p'{t)) + c*(m - 1)1)2(1 + 2a)'^~'^q{t) + eCD 
< e-^/'{-£2N - 1 + c*(m - 1)D^{1 + 2ar~^a) 

+e^\£2N -£L + £c*{m - 1)1)2(1 + 2ay^aL) + £2NK + eCD . 

In view of ()5.12p we get 




(5.18) 



< a' 



1 < (m-l)(l + 2(7)"-2g(i)_ 



(5.19) 



Et<£[ e^\2N - -L) + 2NK + CD] < , 




if L > is sufficiently large. This implies that Ei 
Now, assume that d^{t,x) < so that 



aU'El > 0. 



d'^it, x) — £p{t) 



<-K <-l. 



(5.20) 



e 
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If El < the conclusion Ei > follows. Let us now assume E^ > 0. The 
above study for the case < d^{t,x) < ep{t) implies a fortiori that 

e{Np{t) - p'{t)) + c*(m - l)L>2(i + 2a)"'~^q{t) + eCD < . 

Therefore 

\E*i\<N\d%t,x)-epit)\, 
and we deduce from (|5.20p and (|5.4p that 

1^1 1 < eCiaU. 

Summarizing, we obtain that, in any cases, 

eC'[ut] > [/ {-a + o"^ + [/ (a^ - a"') - eC^a''' + eq'} , 

since a = 1 + g > 1 and with C4 := Ci + C2 + 2C3. Since U < 1, a > 1 and 
m > 2, the inequality — a + a™ + C/ (a^ — a™) > —a + = q + q"^ > q holds. 
Therefore, using |a| < 1 + 2a and substituting the expressions for q{t) and 
q'{t), we see that 

sC^[u+] > U [eaLe^* - eCi{l + 2(7)'" + ae^L'^e^^] 

> Ue{aL-C4{l + 2a)"'} 

>0, 

if L > is sufficiently large. 

This completes the proof of Lemma 15.21 □ 

6 Proof of Theorem 12.2 

By fitting the sub- and super-solutions for the generation into the ones for 
the motion, we are now in the position to prove our main result. 

Let r] G (0, 1/2) be arbitrary. Choose a that satisfies ()5.12p and 

a<|. (6.1) 

By Theorem there exists Mq > such that (gj]), (jM]) and ([O]) hold 
with the constant rj replaced by a/2. Recall that are the sub- and super- 
solutions constructed in (jS.lip . 

Lemma 6.1 (Ordering initial data). There exists K > such that for all 
K > K , all L > 0, all £ > small enough, we have 

u^{0,x) <u%f,x) <uf{0,x) , (6.2) 

for all X G O. 
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Proof. We first prove 

n-(0, x) = {l-a- eaL)U f^^i^i^l±^ \ < u%t',x) . (6.3) 



If X is such that d^(0,x) > —Ke, this is obvious since the definition of U 
then implies u~{0,x) = 0. Next assume that x is such that d'^{0,x) < —Ke. 
Let us denote by d{t,x) the signed distance function to r^. Note that, in 
view of hypothesis (jl.ip . the mean value theorem provides the existence of 
a constant Kq > such that 

if d{0,y)<-Koe then uo{y) > MqE . (6.4) 

Moreover in view of the definition of p^''^''^^* in (|2.4p and the compactness of 
To, there exists Kq > such that, for e > small enough, 

if d''{0,x) < -Koe then d{0,'^{0,f,x))<-Ko£. (6.5) 

Hence, if we choose K > Kq, we deduce from (|6.5|) . (|6.4p and (|4.2p (with rj 
replaced by a/2) that u'^{t^,x) > 1 — f . Since U <1, this proves (|6.3p . 
Next we prove 

u'{f, x)<{l + a + eaL)U J^^li^^i^Jl:^^ = ^+(o, x) . (6.6) 

In view of the definition of p^''^^''^* in (|2.4p and the compactness of Tq, there 
exists Ki > such that 

if d^(0,x) > i^ie then dist($(0, ^^ rr), Oq) > MqE . (6.7) 

If X is such that d'^{0,x) > Kie then it follows from ()6.7p and Theorem 14.11 
{in) that u^(t^,x) = 0, which proves (j6.6p . Next assume that x is such that 
(i^(0,x) < Ki£. Since U is non increasing we have 



f d%0,x) - K£\ . _ a 



(1 + a + eaL)C/ ^ ' ^ J > (1 + a)?7(ii:i - ET) > 1 + - , 

if K > i^i (recall that U{-oo) = 1). Then ([621) follows from (fiT]) (with rj 
replaced by cj/2). □ 

We now prove Theorem 12. 2[ 

Proof. We fix > 1 large enough so that Lemma [6 . 1 1 holds . and L > large 
enough so that Lemma [5.21 holds . Therefore, from the comparison principle, 
we deduce that 

uj{t,x) <u''{t + f,x) <uf{t,x) forO<t<T-f. (6.8) 
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Note that, since 

r 1 iid%t,x) < 
hmuf{t,x) = < (6.9) 

[ ifd%t,x) > 0, 

for t > 0, (|6.8|) is enough to prove the convergence result, namely Corollary 
12.31 We now choose C large enough so that 

{l-^r])U{-C + e^^ + K)>l-r]. (6.10) 
Note that this choice forces 

C>e^'^ + K. (6.11) 

In the sequel we prove (|2.6p . 

Obviously, if e > is small enough, the constant map z"*" = 1 + is a 
super-solution. Therefore we deduce from Theorem l4.1l (i) that u^{t+t^,x) G 
[0,1+7]], for alio <t<T-t^. 

Next we take x € ri^''^"^* \J\fcs(rf^'^^^^), i.e. 

d''{t,x) < -Ce. (6.12) 
For £ > small enough, we have 

uj{t,x) >{l-a-£aLe^^)U{-C + e^^ + K) 

> {l-^a)U{-C + e^^ + K) 

> {l-^rj)Ui-C + e^^ + K) 

where we have successively used ()6.ip and (|6.10p . In view of (|6.8p this 
implies that n^(t + t^,x) >1- r], for ah < t < T - t^. 

Finally we take xe{n\ ftf^''^^) \Afce{Tf'^''^'), i.e. 

d%t,x)>Ce. (6.13) 

Using ()6.1ip we see that, for e > small enough, d^{t,x) — £p{t) > so 
that uf{t,x) = 0, which, in view of ()6.8p implies that u'^{t + t'^,x) = 0, for 
0<t<T-t^. 

This completes the proof of Theorem 12.21 □ 

Acknowledgement. The first author is grateful to Hiroshi Matano for a 
fruitful discussion in Dresden. 
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